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Abstract. The time behaviour of the radial abundance gradients in the galactic disk is investigated on the basis 
of four different samples of planetary nebulae, comprising both smaller, homogeneous sets of data, and larger, 
albeit non-homogeneous samples. Four different chemical elements are considered, namely, oxygen, sulphur, argon 
and neon. Our analysis support our earlier conclusions that, on the average, the radial abundance gradients have 
flattened out in the last 6 to 8 Gyr. 
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1. Introduction 

One of the most interesting observational properties that 
can be devised in order to constrain chemical evolution 
models is the magnitude as well as the space and time 
variations of the radial abundance gradients in the Galaxy 
(see for example Maciel & Costa 120031 Maciel 120001 and 
Henry & Worthey 119991 for recent reviews). Regarding 
the time variation of the gradients, which is possibly the 
most important information that can be obtained from 
the abundance variations in the galactic disk, the galactic 
subsystem of planetary nebulae (PN) plays a particularly 
important role. These objects have relatively well deter- 
mined chemical compositions, which are essentially based 
on the analysis of bright emission lines, and are originated 
from stars within a reasonably large mass - and therefore 
age - bracket while on the main sequence. As a conclu- 
sion, PN can be used to investigate the possibility of a 
time variation of the average gradient in the galactic disk, 
in the same way as the open cluster stars (see for example 
Friel 1995, Friel et al. 2002). 

In the first paper of this series, Maciel et al. H2003I 
hereafter referred to as Paper I) argued in favour of a 
time flattening of the 0/H gradient from roughly —0.11 



Send offprint requests to: W. J. Maciel 

* Based on observations made at the European Southern 
Observatory (Chile) and Laboratorio Nacional de Astrofi'sica 
(Brazil) 



dex/kpc to —0.06 dex/kpc during the last 9 Gyr, or from 
—0.08 dex/kpc to —0.06 dex/kpc in the last 5 Gyr, on the 
basis of a large sample of PN in the galactic disk. More 
recently, Maciel et al. I|2005al Paper II) extended the orig- 
inal discussion by (i) including S/H data, (ii) adopting 
[Fe/H] X 0/H and [Fe/H] x S/H conversions so that an 
estimate of the corresponding [Fe/H] gradient could be 
made, (iii) estimating the average gradient from Cepheid 
data and (iv) taking into account some recent determi- 
nations of the gradients from young objects, such as HII 
regions and stars in OB associations, and presenting a 
general comparison of the derived gradients of all objects 
considered. It has been shown in Paper II that the derived 
results lead to a consistent interpretation of a time flat- 
tening of the radial abundance gradients in the galactic 
disk, in agreement with the main results of Paper I. 

In the present paper, we extend our previous results 
by taking two additional steps: (i) we consider four chem- 
ical elements in the discussion, namely oxygen, sulphur, 
argon and neon, so that the ratios 0/H, S/H, Ar/H and 
Ne/H can be studied; (ii) we take into account four dif- 
ferent PN samples, including three additional sets of data 
apart from the original abundances considered in Papers I 
and II. These samples are generally homogeneous sets, in 
the sense that all objects have been analyzed by the same 
group, under similar conditions, which includes data re- 
duction and abundance determination. This is particularly 
important in order to make the previous results more ro- 
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Table 1. Number of PN in each sample. 



Sample 


O/H 


S/H 


Ar/H 


Ne/H 


Basic Sample 


234 


117 


114 


118 


Henry 


81 


77 


78 


78 


Perinotto 


115 


112 


107 


92 


lAG/USP 


71 


58 


62 


18 



bust, since all the large samples presently available are 
not homogeneous, being compiled from a few individual 
homogeneous samples. As we will show in the next few 
sections, the general conclusions of Papers I and II are 
maintained. 

2. The samples 

2.1. The basic sample 

The basic sample is essentially the same used in Papers I 
and II (Maciel et al. i2003w.2005a|D . and has been described 
in detail in Paper I, to which the reader is referred. It is 
based on an earlier sample by Maciel & Koppen (1994') 
and Maciel & Quireza 1,1999,) . with some additional ob- 
jects presented in Costa et al. H2l)U4|l . This sample cov- 
ers a large range of galactocentric distances, roughly from 
i? = 4 to i? = 14 kpc, with a stronger concentration of 
objects located up to i? = 12 kpc. The main additions to 
our earlier sample are the anticentre objects from Costa 
et al. I|2(J(J4|I . which are typically located at i? > 10 kpc. 
Such a range is clearly wide enough to derive accurate 
gradients, and in fact constitutes one of the strong points 
of using PN to study abundance variations in the galactic 
disk. 

The total number of planetary nebulae for which reli- 
able abundances of O, S, Ar and Ne have been obtained 
is shown in Tabled for the basic sample, as well as for the 
remaining samples considered in this paper. These addi- 
tional samples have approximately the same total galac- 
tocentric distance range as the basic sample, but they are 
smaller and more dispersed, with comparatively less ob- 
jects at i? > 10 kpc. 

It can be seen that the basic sample is the largest sam- 
ple considered, but it should be kept in mind that it is a 
compilation, albeit careful, of several different determina- 
tions in the literature. Therefore, it lacks the homogeneity 
of smaller samples available, so that it is interesting to in- 
vestigate the extent to which the non-homogeneity degree 
affects the results of our analysis. Unfortunately, no sam- 
ple is presently available that includes at the same time a 
high degree of homogeneity and a size large enough to be 
considered a complete sample - e.g. by including objects 
in the whole range of expected abundances. 

The distances of the PN in the basic sample come from 
Maciel & Quireza ifTM^ and Maciel & Koppen ifT^Mjl . 
and are basically from the statistical distance scale by 
Maciel (^1984() . For those nebulae not present in the cata- 



logue of Maciel H1984|l . or for which only distance limits 
were available, we have adopted values from more recent 
scales, especially Cahn et al. H1992|l (see the online mate- 
rial associated with Maciel & Quireza, 119991 and Table 1 
of Costa et al., I2()()4|l . The same sources were used for 
the remaining samples. The problem of the distances to 
PN and their effect on the abundance gradients has been 
discussed elsewhere (see for example Maciel & Koppen, 
^Sf). The analysis of large PN samples is necessarily 
based on statistical distances, which include a much larger 
number of objects than individually determined distances. 
Moreover, the uncertainties of the individual distances are 
often larger than generally assumed, and in many cases 
not very different from those associated with the statisti- 
cal distances. We may conclude that the use of different 
statistical distances may have a large influence on the dis- 
tance of a particular object, but the global effect on a 
large sample is negligible, that is, the distance distribu- 
tion is essentially unaffected. In other words, the abun- 
dance gradients are not sensibly affected, as well as their 
derived time variation. A recent study by Maciel & Lago 
(|2005J led to a similar conclusion regarding PN kinemat- 
ics, which depends on the PN positions - hence distances 
- on the galactic disk. It was found that the use of four 
different statistical distance scales fMaciel [19841 Cahn et 
aLCnni van de Steene & Ziilstra lBMl and ZhanEf [T?M)l 
for disk PN leads to a similar rotation curve, which is not 
very different from the Population I curve derived from 
HII regions fClemens ll985|l . 

2.2. The sample by Henry et al. (2004) 

The first homogeneous sample comes from the work of 
Henry and co-workers (Henry et al. 12004(1 . who derived 
abundances of a set of 85 planetary nebulae based on ob- 
servations at KPNO and CTIO. The detailed discussion of 
each step of the derivation is given in a series of previous 
papers (Kwitter & Henrv l^UHTl Milingo et al HUU^ .h. and 
Kwitter et al. '2003). The abundances have been derived in 
a consistent and homogeneous manner and, especially, the 
sulphur abundances are expected to be particularly well 
determined, as they have included near-IR lines of [SHI] 
at A9069, 9532A, apart from the optical lines of [SII] at 
A6716, 673lA and [SIII] at A6312A. 

Some of the objects in the Henry et al. analysis are 
halo nebulae, and have been removed from our sample, 
as they are not adequate to study disk properties such as 
abundance gradients. 

Henry et al. 1)2004(1 have also analyzed the existence of 
radial abundance gradients from their sample. Although 
no attempt has been made to investigate the time varia- 
tion of the gradients, they were successful in deriving aver- 
age gradients in the range —0.04 to —0.05 dex/kpc for all 
studied elementes, which are O/H, Ne/H, S/H, Cl/H and 
Ar/H. Since their abundances - and to a large extent their 
distances - are completely independent from the sample 
considered in our previous papers, the fact that measur- 
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able and similar gradients have been detected is a nice 
confirmation of the usefulness of disk planetary nebulae 
as a chemical evolution tool. 

2.3. The sample by Perinotto et al. (2004) 

Recently, Perinotto et al. (2004) have presented a detailed 
analysis of the chemical composition of a large sample of 
planetary nebulae, for which they have used the same ex- 
tinction corrections, electron temperature, densities and 
ionic abundance determinations, as well as ionization cor- 
rection factors, including atomic constants and nebular 
models. Since the measured fluxes originated from differ- 
ent sources, this is not strictly an entirely homogeneous 
sample, but its degree of homogeneity is clearly much 
higher than in the case of a simple compilation, in view 
of rigorous criteria adopted for the inclusion of a given 
object and the same procedure applied to the measured 
fluxes. In fact, Perinotto et al. ()2()[)4|l consider their sam- 
ple as a highly homogeneous data set on galactic PN, de- 
termined with realistic uncertainties. Comparison of their 
abundances with some well determined data in the litera- 
ture, such as in the work by Kingsburgh & Barlow (1994) 
generally produce a very good agreement. This sample is 
considerably larger, amounting to over 130 nebulae, some 
of which have been removed from our analysis since they 
may also belong to the halo or have no accurate distances. 

2.4. The I AG/ US P sample 

The last sample considered in this paper comes from our 
own data, which we will call the lAG/USP sample. This is 
a highly homogeneous sample, although relatively small, 
reaching about 70 nebulae, after removing bulge objects 
and PN in the Magellanic Clouds. The observations have 
been made either at the 1.6 m LNA telescope in Brazil 
or with the 1.52 m ESO telescope in Chile. In both cases. 
Boiler and Chivens Cassegrain spectrographs have been 
used. Details on the observations and reduction proce- 
dures can be found in the original papers (see for example 
Costa et al. 120041 and references therein). 

2.5. Comparison of PN abundances 

A comparison of the abundances of the basic sample and 
the additional samples is shown in Figs.HtoElfor oxygen, 
sulphur, argon and neon, respectively. Regarding the sam- 
ple by Henry et al. i2004p . (filled circles) it can be seen 
that the oxygen, argon, and neon abundances generally 
show a very good agreement with the basic sample, while 
part of the sulphur abundances shows some underabun- 
dance relative to the basic sample. In fact, the lower sul- 
phur abundances have been noticed by Henry et al. H2004|l 
when making comparisons with HII region data, and was 
referred to as the "sulphur anomaly" . Several reasons have 
been suggested to explain such a discrepancy (see the dis- 
cussion by Henry et al.), to which we could add the fact 
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bg(0/H)+12 (Basic Sample) 

Fig. 1. A comparison of the 0/H abundances from the 
samples by Henry et al. (2004, filled circles), Perinotto et 
al. (2004, empty circles) and the lAG/USP sample (stars) 
with data of the basic sample. 

that their [SIII] electron temperatures are systematically 
higher than the corresponding [OIII] temperatures, which 
tends to decrease the derived sulphur abundances. A few 
outliers can be observed in Figs. to 01 and it should be 
noticed that some of them, such as Hb 12 and IC418 have 
also been identified as outliers in the plots presented by 
Henry et al. (2004) . 

The sample by Perinotto et al. H2004I empty circles) 
also shows a generally good agreement with the data in 
the basic sample, especially for oxygen (Fig. ^ and neon 
(Fig. 01). Their sulphur abundances, and to a lesser extent, 
argon, are generally lower than in the basic sample (cf. 
Figs. 121 and O, and show the same sulphur anomaly as 
discussed by Henry et al. 1)200411 . 

The smaller lAG sample (stars) shows several nebulae 
directly on the straight line of Figs. Q to 01 which reflects 
the fact that we have adopted several of the results of this 
sample in our merged basic sample. An exception is neon, 
for which our own data is very limited, as can be seen from 
Fig. H or Tabled 

3. The method 

A detailed discussion of the method used in order to esti- 
mate the abundance gradients at different epochs is given 
in Papers I and II (Maciel et al. IMini I2()()5a|l to which 
the reader is referred for details. Basically, the abundance 
gradients in the form (ilog(X/H)/(ii? (dex/kpc) were de- 
termined for X = O, S, Ar and Ne assuming a linear vari- 
ation of the abundances with the galactocentric distance 
R, namely, neglecting possible space variations within 
the galactic disk of the gradients themselves. A value of 
i?o — 8.0 kpc was adopted for the position of the LSR, as 
in the previous papers. 

The samples were divided into age groups, and the ages 
have been estimated in the following way. First, the heavy 
element abundances, that is, 0/H, S/H, Ar/H and Ne/H, 
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Fig. 2. The same as Fig.Hfor S/H. 
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Fig. 3. The same as Fig.[T]for Ar/H. 



. HENRY 






o PERINOTO 






★ lAG 
















o • ° J-"^ 






• « 


• 


• o 

• • 

o o , 


* ° ' 














o 




• 


• 


• 


o 

9 


7.5 




8 8.5 c 




log(Ne/H) 


+ 12 (Basic Sample) 



Fig. 4. The same as Fig. [T] for Ne/H. 

have been converted into [Fc/H] mctalhcities using a well 
determined [0/H] x [Fe/H] relationship for the galactic 
disc, as well as correlations of S/H, Ar/H and Ne/H with 
the oxygen abundance 0/H. These correlations have been 
studied by several people (see for example Henry et al. 
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AGE (Gyr) 

Fig. 5. The adopted age-metallicity relation from 
Edvardsson et al. H1993|l which depends on the galactocen- 
tric distance R. The AMR are plotted for O /H abundances 
and galactocentric distances in the range 4 < i?(kpc) < 
14. The dots show the PN of the basic sample. 



I2UU4|I . and the slopes are very close to unity, reflecting 
the fact that, as a first approximation, all four element 
ratios are good tracers of the interstellar abundances at 
the time of formation of the PN progenitors. From the 
[Fe/H] metallicities, the ages have been determined using 
an age-metallicity relationship (AMR) by Edvardsson et 
al. I)1993|l . which also depends on the galactocentric dis- 
tance. Such a relationship constitutes an improvement rel- 
ative to the solar neighbourhood AMR generally adopted, 
although it should be considered as an approximation as 
applied to the galactic disk. As an example. Fig. shows 
the adopted age-metallicity relation as applied to the O /H 
abundances of the PN in the basic sample. It can be seen 
that the whole age bracket extends from 1 Gyr to 10 Gyr, 
approximately, with a larger concentration of objects with 
ages around 4 to 5 Gyr. 

The adoption of an age-metallicity relation is proba- 
bly our main uncertainty source, even allowing for some 
galactocentric dependence. Although most average AMR 
in the literature are similar to each other (see for instance 
Rocha-Pinto et al.. I2()()()1l . there is still an ongoing dis- 
cussion on their intrinsic dispersion at any given age. An 
alternative approach involves the direct determination of 
the central star masses, and hence ages, from the neb- 
ular chemical composition, as we have done in Paper I, 
with similar results for the time variation of the abun- 
dance gradients. Such a procedure can also be applied to 
the predictions of theoretical models of intermediate mass 
stars. As an example, we have recently considered mod- 
els by Marigo l)2()()l|l and our preliminary results (Maciel, 
Costa & Lago, in preparation) largely support the results 
based on the age-metallicity relation. 

Once the individual ages have been determined, the 
nebulae in each sample have been divided into two age 
groups, one called the younger Group I, and the other the 
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older Group II. As in Paper II, if tj is the upper age limit 
of Group I, all PN progenitors having ages t < tj will 
belong to Group I, while those with ages t > tj belong 
to Group II. Several values of ti have been considered 
in the range 3.0 < i/(Gyr) < 6.0, and for each of these 
values we have calculated the gradients of Groups I and II 
and the corresponding correlation coefficients. Therefore, 
our conclusions do not depend on the adopted age limit 

that is, they are independent of the detailed manner 
according to which the groups are defined. Of course, for 
values of the age limit very close to the limiting values 
of the adopted range, namely i/ ~ 3 or 6 Gyr, either 
Group I or Group II becomes underpopulated, so that the 
corresponding results are less reliable (see Paper II for a 
detailed discussion). 

It should be stressed that the main goal of this project 
is to establish the time variation of the abundance gradi- 
ents, that is, to determine whether the gradients steepen 
or flatten out with time. Therefore, we are not particularly 
concerned with an accurate determination of the absolute 
magnitude of the gradients at any given epoch. 



4. Results and discussion 

4.1. 0/H 

The main results on the time variation of the abundance 
gradients of the O /H ratio are shown in Fig. where we 
plot the derived gradients as a function of the adopted 
age limit tj. These are the best results of all four ele- 
ments, as the average 0/H abundances are higher, better 
determined and the samples are larger. From Fig.|B| it can 
be concluded that the gradients of the younger Group I 
(empty circles connected by lines) are systematically flat- 
ter than the corresponding gradients of the older Group II 
(filled circles connected by lines), supporting our earlier 
conclusions of Papers I and II. These results apply to all 
four samples considered, especially the larger basic sam- 
ple, the sample by Henry et al. l|^nM| and the lAG/USP 
sample. For the basic sample (top left panel) the results 
are essentially the same as presented in Fig. 1 of Paper II. 

In the case of the sample by Perinotto et al. (fSDM', 
bottom left panel), the gradients of both groups are similar 
for age limits in the range 3 < t/(Gyr) < 4.5, but for larger 
age limits the separation of the two groups is clear. 

In view of the comment at the end of Sec. 3, for 
~ 3 Gyr the younger Group I is underpopulated, and 
the same occurs at t/ ~ 6 Gyr for Group II. Therefore, 
near these limiting values the statistical results are less 
reliable, as the space distribution of the nebulae may not 
be homogeneous. We would then expect the best results 
to be associated with groups containing a reasonably large 
fraction of the whole sample. As an illustration, the age 
limit at which both groups have the same number of el- 
ements (that is, half of the sample) occurs typically for 
3.9 < t/(Gyr) < 4.6, showing no systematic differences 
between the space distributions of the PN groups. 



A detailed discussion of the errors involved in the de- 
termination of the gradients, as well as a statistical anal- 
ysis of the different object samples considered have been 
recently given by Maciel et al. I|2005b|l . For the sake of 
completeness, we have included at the lower right corner 
of the last panel of Fig. an average error bar concerning 
the 0/H gradient from PN. The same procedure has been 
adopted in the corresponding figures of the other element 
ratios, namely S/H, Ar/H and Ne/H. 

Another important quantity that can be used in the 
interpretation of Fig. El is the correlation coefficient of the 
linear fits. For 0/H, the coefficients of all fits shown in 
Fig.iniare large, typically in the range 0.50 < \r\ < 0.90, so 
that the flattening of the gradients of the younger Group I 
relative to the older Group II is real. 



4.2. S/H 

For the S/H ratio, the same behaviour of the 0/H gradi- 
ents is observed, as can be seen from Fig. |7| This is par- 
ticularly true for the basic sample, the sample by Henry 
et al. I|2()04(l and the lAG/USP sample. Again the gradi- 
ents of the younger Group I are systematically flatter than 
for the older Group II, and the correlation coefficients are 
also generally large, in the range 0.50 < \r\ < 0.80. For 
the basic sample, in particular, the correlation coefficients 
of both groups are large, 0.70 < |r| < 0.80. The results of 
the top left panel are also essentially the same as in Fig. 2 
of Paper II. For these three samples, the average gradi- 
ents of Group I are essentially independent of the adopted 
age limit, while the corresponding gradients of the older 
groups are increasingly steep, although the derived values 
for age limits close to the upper limit of 6 Gyr are less 
reliable, as mentioned. 

It is interesting to notice that the behaviour of S/H 
gradients of the basic sample and the Henry et al. sample 
are very similar, as shown by the top panels of Fig. [7| As 
mentioned in Sec. 2, the sulphur abundances of the sample 
by Henry et al. are expected to be very accurate, providing 
a further confirmation of the flattening of the gradients 
for the younger Group I. It can also be concluded that 
the "sulphur anomaly" present in the Henry et al. sample 
does not affect the picture of the time evolution of the 
gradients as derived from the basic sample. 

In the sample by Perinotto et al. H2004|l . S/H shows 
a similar behaviour as 0/H, in the sense that the gra- 
dients of both groups are of the same order of magni- 
tude for smaller values of the age limit, tj < 4.5 Gyr, 
presenting some flattening of the younger Group I at 
higher values of tj. For this sample, the correlation co- 
efficients reach the lowest values, typically \r\ ~ 0.4 to 
0.5. Again for illustration purposes, in all cases of Fig. |7| 
both Groups I and II are approximately of the same size 
for 4.1 < ^/(Gyr) < 4.6. 



6 



Maciel et al.: Time variation of radial abundance gradients III. 



d 



LO 

CN 

O 

^ d 

BASIC SAMPLE ° 
II 


HENRY 

II 


3 4 5 6 

LO 

in 

PERINOTTO ° 
II 


3 4 5 6 

I 

lAO/USP 

II 



3 4 5 6 3 4 5 



Age limit (Gyr) 

Fig. 6. Time variation of the 0/H gradient from planetary nebulae. The four PN samples were divided into two age 
groups, Group I ("younger"), with ages lower than the age Hmit tj, and Group II ("older"), with ages higher than 
tj. The plot shows the 0/H gradient (dex/kpc) of each group as a function of the upper age limit of Group I, ti, for 
the four samples considered. The gradients of the younger Group I (empty circles connected by lines) are flatter than 
those of the older Group II (filled circles connected by lines). An average error bar is given at the lower right corner 
of the last panel. 



4.3. Ar/H 

The Ar/H gradients shown in Fig. |S1 also indicate some 
flattening of the younger Group I in the basic sample, 
especially for age limits tj > 4.5 Gyr. The corresponding 
correlation coefficients are moderate, typically |r| ~ 0.50. 

For the lAG/USP sample, the flattening is more evi- 
dent, as can be seen from the bottom right panel of Fig.|Hl 
where Group I (empty circles) indicate less pronounced 
gradientes than Group II (fllled circles). Here the correla- 
tion coefficients are higher, in the range 0.50 < |r| < 0.80. 

In the sample by Henry et al., some flattening of 
Group I can also be noticed for tj > 4 Gyr. In all pan- 
els of Fig. El both groups have similar sizes for 3.8 < 
t/(Gyr) < 4.6, so that the observed flattening is proba- 
bly real. However, the correlation coefficients in this range 
are modest, typically |r| ~ 0.40 for Group I and somewhat 
higher (|r| ~ 0.60) for Group II, so that this conclusion 
cannot be stressed too much. 

For this element, the sample by Perinotto et al. does 
not lead to any definite conclusions. In this case, the gra- 
dients of both groups are very similar, within the expected 
uncertainties of the gradients themselves, which are typi- 
cally 0.02 dex/kpc, as can be seen from Papers I and II. 
Since the correlation coefficients are modest (|r| ~ 0.4 to 



0.6), we should probably conclude that no differences can 
be observed between the gradients of the different groups. 

4.4. Ne/H 

Finally, Fig.Elshows the results of the Ne/H ratio for our 
adopted samples. We can see that a well defined flattening 
of the younger Group I relative to the older Group II is 
apparent for the Henry et al. H2004|l sample, for which 
the correlation coeflicients are also moderate to large, 
being in the range 0.40 < \r\ < 0.60 for Group I and 
0.50 < |r| < 0.80 for Group II. Such a result is partic- 
ularly important, as the remaining two samples do not 
allow clear conclusions to be drawn. For this sample, both 
groups have approximately the same size for ti ~ 4.3 Gyr. 

The Perinotto et al. sample displays essentially sim- 
ilar gradients for both age groups. The correlation co- 
efficients are large in this case, typically in the range 
0.50 < |r| < 0.80, so that we conclude from this sample 
that the gradients of both groups are similar. 

In the case of the basic sample, acceptable gradients 
can be obtained only for the younger Group I, for which 
the correlation coefficients are well defined, and typically 
in the range 0.40 < |r| < 0.80, so that the derived gradi- 
ents are probably correct. However, for the older groups 
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Fig. 7. The same as Fig.Elfor S/H. 



the gradient is very flat for most of the range in age hm- 
its, and the correlation coefficients are very low, typically 
\r\ < 0.20, which is consistent with the very low gradients 
displayed by this group. This means that the gradients are 
not well defined in this group, so that no comparison of 
the time evolution of the gradients can be made. 

As mentioned before, the lAG/USP sample includes 
only a very limited number of PN with well determined 
Ne/H abundances, so that no reliable gradients can be 
obtained from this sample. The corresponding panel in 
Fig.inishows some indication of a flattening for the younger 
Group I, but the correlations are rather poor and the sam- 
ple is very small. The gradients of the older Group II are 
better, so that we can conclude the results of this sam- 
ple are consistent with the results of the sample by Henry 
et al. It seems that for this element the homogeneity of 
the sample is especially important, as the the less homo- 
geneous samples (the basic and the Perinotto samples) 
present uncertain and even conflicting results, while the 
the homogeneous sample by Henry et al. shows 
very clear results. 

5. Time evolution of the [Fe/H] gradient 

Taken together, the results shown in Figs. El to El gener- 
ally confirm the main conclusions of Papers I and II, in 
the sense that there is a clear tendency for the gradients 
to flatten out, at least during the lifetimes of the older 
objects considered here. Adopting 13.6 Gyr for the age of 
the Galaxy, as in the previous papers, the gradients have 



been flattening out in the last 6 to 8 Gyr, approximately, 
or since the Galaxy was approximately 6 Gyr old. From 
the data considered here, as in the previous papers, noth- 
ing can be said about the time evolution of the abundance 
gradients during the early epochs of galactic evolution, 
that is, for t < 6 Gyr. These results cannot be taken as 
definitive, in view of the uncertainties and assumptions 
considered, especially regarding the abundances and age 
determinations. However, as we will see below, we have 
considerably expanded our original analysis to include new 
elements, samples and galactic objects, so that our con- 
clusions are now based on a more solid foundation. 

The inclusion of Ar and Ne generally supports the 
earlier conclusions drawn from 0/H (Paper I) and S/H 
(Paper II), although these elements are clearly more lim- 
ited in view of the lower, more uncertain abundances, 
and the smaller samples. In this respect, The homoge- 
neous sample by Henry et al. H2()()4|l is particularly use- 
ful, especially for neon. Adopting this sample as repre- 
sentative of the the Ne/H gradient, we would estimate 
for Group I (ages between and 4 Gyr) an average 
gradient of d log (Ne/H) /di? ~ — 0.06 dex/kpc, and for 
Group II (ages between 4 and 8 Gyr) an average gradient 
of c? log(Ne/H) /c?_R ~ —0.11 dex/kpc, which would corre- 
spond approximately to an [Fe/H] gradient d[Fe/ll]/dR ~ 
-0.07dex/kpc and d[Fe/H]/di? ~ -0.13dex/kpc for 
Groups I and II, respectively (see Paper II for details in 
the conversion of observed abundance gradients to [Fe/H] 
gradients). 



8 



Maciel et al.: Time variation of radial abundance gradients III. 




BASIC SAMPLE 

\ L 

3 4 




PERINOTTO 



Fig. 8. The same as Fig.Elfor Ar/H. 



o - 
o 




HENRY 




lAG/USP 



3 



Age limit (Gyr) 



For Ar/H, using the same procedure, and taking 
into account the basic sample as well as the sample by 
Henry et al. H2UU4fl . we would obtain (ilog(Ar/H)/(ii? ~ 
-0.04dex/kpc or d[Fe/H]/di? ~ -0.05dex/kpc for 
Group I and d log ( Ar/H) /di? ~ -0.075 dex/kpc or 
d[Fe/H]/di? ~ -0.09 dex/kpc for Group H. 

These estimates are plotted in FigE]| which repro- 
duces Fig. 8 of Paper H with the following additions: (i) 
the Ar/H gradients from PN converted into [Fe/H] gra- 
dients are included, as described above (empty squares); 
(ii) the same for Ne/H gradients (open triangles); (iii) re- 
cent results by Esteban et al. I|2005|l for 0/H gradients in 
HH regions based on abundances derived from recombina- 
tion lines, converted into [Fe/H] gradients as in Paper H 
(star of David) , and (iv) predictions of theoretical models 
(model A) by Chiappini et al. H2001|l are shown as a dot- 
ted line. The remaining results plotted in Fig. ^1 are the 
same as in Fig. 8 of Paper H, namely: (i) 0/H data from 
PN converted into [Fe/H] gradients, considering three age 
groups (young, intermediate, old) (filled squares); (ii) the 
same for S/H (filled triangles); (iii) [Fe/H] data from open 
clusters from the sample by Friel et al. I|2()()2I filled cir- 
cles); (iv) the same for the sample by Chen et al. H2()();-{l 
open circles); (v) [Fe/H] data from cepheids (Andrievsky 
et al.Enn3ibc,Ennil Luck et al^MJ^, filled star); (vi) OB 
stars and associations from Daflon & Cunha H2004L thick 
cross), and (vii) HH region data from the recent litera- 
ture, as discussed in Paper II (x sign). Finally, the dashed 
line shows predictions by theoretical models by Hou et al. 



H2000(l . which can be compared with the previously men- 
tioned results by Chiappini et al. H2()()l|l as illustrations 
of theoretical models of chemical evolution. The dotted 
vertical line indicates the adopted age of the galactic disk, 
ta = 13.6 Gyr. 

As mentioned in Paper II, our purpose here is to pro- 
vide observational constraints to chemical evolution mod- 
els, rather than to discuss in detail the different models 
available in the literature. Nevertheless, it is interesting to 
comment on the large discrepancy of the models by Hou 
et al. I)2000|l and Chiappini et al. (2Q01J as observed in 
Fig. especially at earlier epochs of the galactic life- 
time. Hou et al. l|2()()()|l adopt an exponentially decreasing 
infall rate and an "inside-out" scheme for the formation 
of the galactic disk, in which a rapid increase of the metal 
abundance at early times in the inner disk leads to a steep 
gradient. As the star formation migrates to the outer disk, 
metal abundances are enhanced in that region, so that 
the gradients flatten out. Similar results have also been 
obtained by Alibes et al. (|2001fl . among others. In the 
models by Chiappini et al. H2001(l . two infall episodes are 
assumed to form the halo and the disk. The latter also has 
an "inside-out" formation scenario, in which the timescale 
is a linear function of the galactocentric distance. The de- 
tailed time variation of the gradients depends on the par- 
ticular model, but in general some steepening of the gra- 
dients is predicted, which docs not seem to be supported 
by the data presented in Fig. EH Possibly, the different 
timescales for star formation and infall in the models by 
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Hou et al. (^lOT) and Chiappini et al. fWU^ are respon- 
sible for their different predictions on the time variation 
of the abundance gradients. 

It can also be seen that the new results fit nicely the 
picture previously sketched in Paper II, so that the same 
average flattening rate can be advocated here, namely, 
d[Fe/H]/di? - 0.005 to -0.010 dex kpc"! Gyr"! for the 
last 6 to 8 Gyr. This value refers to a spatially- averaged 
gradient, corresponding to a linear variation of the abun- 
dances with the galactocentric distance. In practice, the 
situation is more complicated, and all quoted theoreti- 
cal models predict some space variations of the gradients 
along with their temporal variations. Since the observa- 
tional results on the space variations are still very limited 
and often contradictory, further work is needed in order 
to definitely establish the characteristics of this important 
observational constraint. 

Acknowledgements. This work was partly supported by 
FAPESP (grant 02/08816-5), CNPqand CAPES. Observations 
at ESO/Chile were possible through the FAPESP grant 
98/10138-8. 

References 

Alibes, A., Labay, J., & Canal, R. 2001, A&A, 370, 1103 
Andrievsky, S.M., Bersier, D., Kovtyukh, V.V., Luck, R.E., 

Maciel, W. J., Lepine, J.R.D., & Beletsky, Yu. V. 2002b, 

A&A, 384, 140 

Andrievsky, S.M., Kovtyukh, V.V., Luck, R.E., Lepine, J.R.D., 
Bersier, D., Maciel, W.J., Barbuy, B., Klochkova, V. G., 
Panchuk, V. E., & Karpischek, R. U. 2002a, A&A, 381, 32 



Andrievsky, S.M., Kovtyukh, V.V., Luck, R.E., Lepine, J.R.D., 
Maciel, W.J., & Beletsky, Yu. V. 2002c, A&A, 392, 491 

Andrievsky, S.M., Luck, R.E., Martin, P., & Lepine, J.R.D. 
2004, A&A, 413, 159 

Cahn, J. H., Kaler, J. B., & Stanghellini, L. 1992, A&AS, 94, 
399 

Chen, L., Hou, J. L., & Wang, J. J. 2003, AJ, 125, 1397 
Chiappini, C, Matteucci, F., & Romano, D. 2001, ApJ, 554, 
1044 

Clemens, D. P. 1985, ApJ, 295, 422 

Costa, R. D. D., Uchida, M. M. M., & Maciel, W. J. 2004, 

A&A, 423, 199 
Daflon, S. & Cunha, K. 2004, ApJ, 617, 1115 
Esteban, C, Garcia-Rojas, J., Peimbert, M., Peimbert, A., 

Ruiz, M. T., Rodriguez, M., & Carigi, L. 2005, ApJ, 618, 

L98 

Edvardsson, B., Andersen, J., Gustafsson, B., Lambert, D.L., 
Nissen, P.P. & Tomkin, J. 1993, A&A 275, 101 

Friel, E. D. 1995, An. Rev. A&A, 33, 381 

Friel, E. D., Janes, K. A., Tavarez, M., Scott, J., Katsanis, R., 
Lotz, J., Hong, L., & Miller, N. 2002, AJ, 124, 2693 

Henry, R. B. C, Kwitter, K. B., & Balick, B. 2004, AJ, 127, 
2284 

Henry, R. B. C. & Worthey G. 1999, PASP, 111, 919 
Hou, J. L., Prantzos, N., & Boissier, S. 2000, A&A, 362, 921 
Kingsburgh, R. L., & Barlow, M. J., 1994, MNRAS, 271, 257 
Kwitter, K. B., & Henry, R. B. C, 2001, ApJ, 562, 804 
Kwitter, K. B., Henry R. B. C, & Milingo, J. B., 2003, PASP, 
115, 80 

Luck, R.E., Gieren, W. P., Andrievsky, S.M., Kovtyukh, V.V., 

Fouque, P., Pont, F., & Kienzle, F. 2003, A&A, 401, 939 
Maciel, W. J. 1984, A&AS, 55, 253 

Maciel, W. J. 2000, in The Evolution of the Milky Way ed. F. 
Matteucci & F. Giovannelli (Dordrecht: Kluwer), 81 



10 



Maciel et al.: Time variation of radial abundance gradients III. 



o 
d 



u 

CL 



a: 



■ Planetary Nebulae (0/H data) 

A Planetary Nebulae (S/H data) 

□ Planetary Nebulae (Ar/H data) 

A Planetary Nebulae (Ne/H data) 



• Open Clusters (Friel) 
o Open Clusters (Chen) 
•k Cepheids 

* OB stars 

X HII regions Esteban 



Chiappini 



Hou 




t„ = 13.6 Gyr 



10 



12 



14 



Time (Gyr) 

Fig. 10. Time variation of the [Fc/H] abundance gradient (dex/kpc). The sme as Fig. 8 from Paper II, with the 
following additions: (i) Converted [Fe/H] gradients from Ar/H abundances in FN analyzed in the present paper 
(empty squares); (ii) the same for Ne/H gradients (open triangles); (iii) recent results by Esteban et al. f2005(l for 
HII regions (star of David), and (iv) theoretical tracks from Chiappini et al. (POOl, dotted line). The remaining data 
are: (i) [Fe/H] gradients from PN calculated from 0/H data (filled squares), (ii) the same for S/H (filled triangles); 
(iii) Open cluster data from Friel et al. I|20n2l filled circles); (iv) the same from Chen et al. pOO^J empty circles); (v) 
Cepheid data (filled star); (vi) OB stars in associations (thick cross), (vii) HII regions from the literature (x sign), 
and (viii) theoretical models by Hou et al. 1)20001 dashed line). The dotted vertical line shows the adopted age of the 
galactic disk, Iq — 13.6 Gyr. 
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